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Abstract. This paper considers random (non-Hermitian) circulant matrices, and proves several 
results analogous to recent theorems on non-Hermitian random matrices with independent entries. 
In particular, the limiting spectral distribution of a random circulant matrix is shown to be complex 
normal, and bounds are given for the probability that a circulant sign matrix is singular. 



1. Introduction 

Given a sequence Xo, Xi, . . . of independent complex random variables, denote by C n the n x n 
random circulant matrix with first row Xo, . . . , X n _i: 
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It is well-known (and easy to verify) that the eigenvalues of C ra are 



n-1 



k = 0, 



where uj n = e 2 " Kt l n . The eigenvalues of random circulant matrices have only recently been studied 
explicitly in the literature, e.g. in (UEJIH], which consider real symmetric circulant matrices and 
variants of them. Other models of random matrices lying in some linear subspace of matrix space 
have also been studied recently in [U [71 EJ [HI [19] , among others. These papers mainly prove results 
analogous to classical theorems for Wigner-type random matrices, i.e. symmetric (or Hermitian) 
random matrices whose entries are all independent except for the symmetry constraint. 

This paper mainly considers the eigenvalues of the random circulant matrices S n with no symme- 
try constraint, and prove results analogous to recent theorems about random matrices with all inde- 
pendent entries. In Section [2] we investigate the limiting spectral distribution of large-dimensional 
random circulant matrices. In Section[3]we consider the joint distribution of eigenvalues of circulant 
matrices with Gaussian entries and observe some consequences, especially for the distributions of 
extreme eigenvalues. Finally, in Section 0] we investigate the probability that a circulant matrix 
with ±1 Bernoulli entries is singular. 



2. Limiting spectral distribution 

We denote by 



1 n— 1 

(2.1) A ^ = ^E^. 



j=0 
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for k = 0, . . . , n — 1 the eigenvalues of n 1//2 6 n , and by 



1 n— 1 



fc=o 



the empirical spectral measure of n 1//2 C n - Theorems [U and [2] show that as n — > oo, \i n converges 
to a universal limiting measure, namely the standard complex Gaussian measure 7c with density 
le - ' 2 ' with respect to Lebesgue measure on C. Theorem Q] shows that this convergence holds in 
probability under very weak assumptions on the random variables Xj (including in particular the 
case of i.i.d. random variables with finite variance); Theorem [2] strengthens this to almost sure 
convergence under much stronger assumptions. 

Theorems Q] and [2] are analogues for circulant matrices of the circular law for eigenvalues of 
random matrices with independent entries, which was recently proved by Tao and Vu in [25] in the 
almost sure sense for i.i.d. entries with finite variance. It is likely that the conclusion of Theorem 
[2] also holds in this level of generality. However, the result of [25J and other results leading up to it 
were made possible by recent advances in controlling how close a random matrix with i.i.d. entries 
is to being singular. As discussed in Section |4] below, random circulant matrices are frequently 
much more likely to be singular; thus it may be difficult to prove the optimal result for circulant 
matrices. (On the other hand, the fact that circulant matrices are normal may make such a result 
approachable via the classical moment method used extensively for Hermitian random matrices.) 



Theorem 1. Suppose that the Xj satisfy 
(2.2) 

for some a £ C and 



EXj = 0, EXj =a, E \Xj\ 2 = 1, 



ra-1 



(2.3) 



n— >oo n 



j=0 



for every e > 0. Then \x n converges in expectation and weak-* in probability to j£, the standard 
complex Gaussian measure on C. 

Observe that the hypotheses cover the case of i.i.d. complex random variables with finite variance, 
as well as real random variables and rotationally invariant distributions satisfying the Lindeberg 
condition (12.31). 



Proof. Without loss of generality we may assume a£l, Observe that for a measurable set A C C, 



1 n— 1 -. n—1 

E^) = TF[A t ei] = -V] 



k=0 



k=0 



-. n—1 



E(Re utfXj ) 2 E(Re utfXj ) (Im u^Xjf 



E(ReLo n K X j )(Imu; n K X j 



E(Smu4! l> X j ) 



We consider as a sum of independent random vectors in IR 2 = C, so we will need 
Cov^Xj] 

The identities 
(2.4) 



2(10 + w)z = (Rew)(Re z) + i(Rew)(lmz) 



h(w — w)z = -(Imw)(Imz) + i(Reu))(Imz). 
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for w, z G C are useful. Letting w = z = w„ Xj, 

n—l 1 n— 1 -. 

J>(Re^%-) 2 = -Re^E^X 2 + |A,| 2 ) = - Re J>u^ fc + 1). 



n-l 



i=o j=o i=o 



Since u; 2fc is an nth root of unity, X^j=o = ® unless ui 2k = 1, which is the case only if k = or 
fc = n/2. Thus unless k = or n/2, 

n-l 



i£E(Re^) 2 = i. 



j=0 



The other covariances are computed similarly and it follows that for k ^ 0, n/2 



1 n_1 1 
-£CovK%.) = -/ 2 . 



Since |w n Xj\ = \Xj\, by f|2.3|) we can now apply a quantitative two-dimensional version of 
Lindeberg's central limit theorem to the complex random variables {ub Xj | < j < n — 1}. Let 
iCCbe measurable and convex and assume k ^ 0, n/2. By the proof of Corollary 18.2]), there 
is a function h(n) with lim n _ >o0 h(n) = 0, which depends on A and the rate of convergence in (|2,3p 
but is independent of k, such that 



1 n— 1 

3=0 



Therefore 

(2.5) \EnniA) - 7C (A)| < - |P[A fc 7C (A)| < 



n-l 



< /i(n). 



2 + (n-2)/i(n) n - 



0. 



n 
fc=o 

Thus the claimed convergence holds in expectation 
Next observe that 

M- 

E^ n {A) 2 = E ' 



n 



- n— 1 

^2 W^eA 



n-l 



^ Y, r[(\ k ,x e )eAxA]. 



kl=0 



As before, we can compute the sums of the relevant covariance matrices using ()2.4|) . in this case 
with w = uJnXj and z = LOnXj. We obtain 

n-l 



1 1 

-^CovK%-,^) = -/4 



3=0 

except when k = or n/2, £ = or n/2, k + ^ = or n, or fc = £. These exceptional cases 
account for fewer than 6n of the n 2 possible values of (k,£). Applying the four-dimensional case of 
Lindeberg's theorem similarly to above, we have 

1 n— 1 

(2.6) \Eu n (Af- lc (A) 2 \ l P [( A ^) £AxA] - 7 c(A) 2 | ^> . 



k,l=0 



We are of course using here that standard Gaussian measure on C 2 is the two- fold product of 7c. 
Prom ([23]) and (pUl) it follows that 



E[» n (A) - 1C (A)] 2 = [Eu n (Af - lc (Af] - 2 lc (A)[E^ n (A) - lc (A)] 



0. 
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Thus for every convex measurable A C C, the random variable n n (A) converges to 7c (A) in I? 
and hence in probability. □ 

Bose and Mitra [3j proved earlier a result analogous to Theorem [1] for real symmetric circulant 
matrices, which is a circulant matrix analogue of Wigner's semicircle law; the limiting distribution in 
this case is a real Gaussian measure. This result was strengthened from convergence in probability 
to almost sure convergence by Massey, Miller, and Sinsheimer [19] . The main result of [3] is a 
similar result for circulant Hankel matrices, which amounts to studying the singular values instead 
of eigenvalues of S n (see Remark 1.2] about the relationship between Hankel and Toeplitz 
matrices). Since circulant matrices are normal, their singular values are simply the moduli of their 
eigenvalues; thus this latter result is essentially a corollary of Theorem [1] The proof of Theorem Q] 
follows the basic outline of the proofs of which assume the Xj are i.i.d. with finite third absolute 
moments; the greater generality of Theorem [1] is achieved by applying Lindeberg's theorem, instead 
of the Berry-Esseen theorem as in [2] . 

The statement of Theorem [1] assumes that the same sequence of random variables Xj is used to 
construct C n for every n. The proof shows however that the result generalizes directly to circulant 
matrices constructed from a triangular array of random variables. The same comment applies to 
Theorem [2] below. 



Theorem 2. Suppose that, in addition to (|2.2p one of the following holds. 

(1) There exists a K > such that for each j, \Xj\ < K almost surely. 

(2) There exists a K > such that for each j, the distribution of Xj satisfies a quadratic 
transportation cost inequality with constant K (see below for the meaning of this). 

Then fj, n converges weak-* to 7c almost surely. 

Recall that a probability measure fi on M. d is said to satisfy a quadratic transportation cost 
inequality with constant K > if 



inf / / \x-y\ A dvr(x,y) < ^2KH(^\v) 
■K&n{ti,v)J J 

for every probability measure v on M. d . Here n(//, v) is the class of probability measures on M. d x M rf 
with marginals jx and v respectively, and H{u\v) is relative entropy. Such an inequality is satisfied 
in particular if \i satisfies a logarithmic Sobolev inequality. See [181 Chapter 6] for background and 
references. 

Proof. We begin with the assumption that a quadratic transportation cost inequality is satisfied. 
This assumption implies (and is essentially equivalent to, see [12] ) the following concentration 
inequality for Lipschitz functions of (X±, . . . , X n ). Let F : C n — ► R have Lipschitz constant l-F^p = 
L. Then 

(2.7) P[ \F(X!, ...,X n )- EF(Xi, . ..,X n )\>t]< Ce^ 

for every t > 0, where c, C > are absolute constants. Now for a compactly supported smooth 
function / : C — ► R, define 



n 



1 

F{xi,...,x n ) = ~y2f(xj) 
so that F(X\, . . . , X n ) = f c f dfj, n . Then 

\F(xi,...,x n ) - F(y 1 ,. . . ,y n )\ < ] -^ 2 -^2\x j - yj\ < 



I /Lip 1 1 , I /Lip 
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so | J f r | Lip < n l l 2 I/Iljp- Combining Theorem Q3 ()2.7|) . and the Borel-Cantelli lemma, we obtain 
that for each compactly supported smooth function / : C — > R, 



/ dfj, n ► / / d-yc 

Jc 

almost surely. Applying this for a countable dense family of such / proves the theorem. 

The case of bounded entries is treated similarly using Talagrand's famous concentration inequal- 
ity for convex Lipschitz functions of bounded random variables [23] (also see [18} Section 4.2]); in 
this case an extra step is required to handle non-convex test functions / : C — > R. 

As above, let / : C — > R be compactly supported and smooth. Let R > be such that f(x) = 
for |x| > R and let A > be such that the eigenvalues of the Hessian D 2 f are bounded below by 
—A. For x G C define 



9{x) 



| \x\ 2 if |x| < R, 



[XR(\x\ - f) if \x\ > R. 

Then l^lLip = an( ^ so 1/ + fflLip — 1/lLip Furthermore, since D 2 (f + g) > 0, both g and 

/ + g are convex. Applying the above argument to / + g and g using Talagrand's concentration 
inequality in place of (|2.7p implies that with probability 1, 



/ (/ + 9) dfi n "^°°> / (/ + g) d^c and / g dfi n "^°°> / g dj c 
Jc Jc Jc Jc 

and hence 

/ / d(i n """°°> / / d-yc- 
Jc Jc 

Again, applying this for a countable dense family of such / proves the theorem. □ 



We remark that a slight generalization of the argument in the last paragraph shows that convex 
Lipschitz functions on Mr form a convergence-determining class for the family of probability mea- 
sures on M d with respect to which such functions are integrable. This fact is presumably well-known 
to experts but we could not find a statement of it in the literature. 



3. Gaussian circulant matrices and extreme eigenvalues 

The following result is a circulant analogue of classical formulas (found, e.g., in [21]) for the joint 
distribution of eigenvalues of random matrices with independent complex Gaussian entries. 

Proposition 3. Let each Xj have the standard complex normal distribution. Then the sequence 
Ao, • • • , A n _i of eigenvalues of n~ l l 2 Q n is distributed as n independent standard complex normal 
random variables. 

Proof. The map (Xq, . . . , X^-i) i— > (Ao, . . . , A n _i) defined by (|2.1|) is easily checked to be a unitary 
transformation of C n , so it preserves the standard Gaussian measure on C n . □ 



An easy consequence of Proposition is that in this setting, E/x n = 7c for every n, and not only 
in the limit n — * 00 as guaranteed by Theorem [TJ 

If each Xj is a standard real normal random variable, then the same observation implies that 
the sequence of eigenvalues Aq, . . . , A n _i are jointly Gaussian random variables, but with singular 
covariance since this sequence will lie in an n-dimensional real subspace of the 2n-dimensional 
space C n . On the other hand, Proposition [3] does have the following simple analogue for complex 
Hermitian circulant matrices with Gaussian entries. 
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Corollary 4. Let C„ be a Hermitian random circulant matrix 

' Yp Yi Yn-2 Y n -{ 

Y n -! Y Y 1 y„_ 2 
i Y^[ Y \ 

% _ _Yo_ Fi 

. *l Y 2 y n _x Y _ 

where Yq, . . . , Y\ n i2\ are independent, Yj = Y n —j for j > n/2, Yq and Y n /2 (if n is even) have 
the standard real normal distribution, and Yj has the standard complex normal distribution for 
1 < j < n/2. Then the eigenvalues of n~ l l 2 Qn are distributed as n independent standard real 
normal random variables. 

Proof. If C n is as in Proposition [3J then has the same distribution as -^(S n + C* ). Because C n 
is normal, the eigenvalues of n _1 ' 2 C^ are 

-^(\k + h) = V2Re\ k 

for k = 0, . . . , n— 1, which by Proposition [3] are independent standard real normal random variables. 

□ 



H 



The Gaussian Hermitian circulant matrix of Corollary 0] should be thought of as a circulant 
analogue of the Gaussian Unitary Ensemble, which (up to a choice of normalization) is defined as 
-^(G + G*), where G is an n x n random matrix whose entries are independent standard complex 
normal random variables. 

Corollary 5. Let C n and be as in Proposition and Corollary [^} Let ct\ > • • • > a n > be 

the eigenvalues o/n~ 1 C n C* and let f3\ > ■ ■ ■ > [3 n be the eigenvalues of 'n _1//2 C^ . Then a n has an 
exponential distribution with mean 1/n, and 

ct\ — log n 



and 



Py\ (r ATI loglogn + log4vr\ 

V21ogn Px - V21ogn 

V 2v21ogra / 



both converge in distribution as n — > oo to the Gumbel distribution with cumulative distribution 
function e~ e x . 

Proof. Observe that since S n is normal, the eigenvalues of n _1 C n C* are \Xk\ 2 (though not usually 
in the same order), which by Proposition [3] are distributed as independent exponential random 
variables. The corollary follows by combinging this fact and Corollary [J] with classical theorems of 
extreme value theory (see [T71 Chapter 1]). □ 

The asymptotic distributions of Pi and ct\ are circulant matrix analogues of famous results of 
Tracy and Widom [26] and Johnstone [15], respectively. Davis and Mikosch [TO], who did not 
consider random circulant matrices explicitly, proved (with slight modifications) what amounts to 
a universality result for a±, which shows the conclusion follows for quite general distributions of 
the Xj. Following their method Bryc and Sethuraman [8] proved an explicitly stated universality 
result for these are then the circulant matrix analogues of the results of Soshnikov in [231 [22] 
respectively. The rough orders of magnitude of a\ and fi\ follow in even greater generality from 
work of the author [20] and Adamczak [1] (see the remarks in }20j, Section 3.1]). 
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4. Singularity of circulant sign matrices 

In this section we specialize to the case in which ¥[Xj = — 1] = P[Xj = 1] = 1/2 for every j and 
consider the probability that C n is singular. The corresponding problem for random nxn matrices 
M n with independent ±1 entries has a long history. An old conjecture (see [16]) claims that 

P[M n is singular] = Q + o(l) 

which is asymptotically the probability that two rows of M n are equal up to sign. The best result 
currently known, proved by Bourgain, Vu, and Wood in [6], is 

F[M n is singular] < f-^= + o(l 

By contrast, the following result shows that the singularity probability of an n x n random circulant 
matrix with ±1 entries depends strongly on the number-theoretic properties of n. 

Theorem 6. Let F[Xj = -1] = F[Xj = 1] = 1/2 for each j. If n is even, then 
(4.1) -^L < P[6 n is singular] < 



where c\,c<i > are absolute constants. If n > 3 is odd, then 
(4.2) P[C« is singular] < min < C3 



d(n) \ - 




n 

where C3 > is an absolute constant, d(n) is the number of divisors of n, and <p(m) is the number 
of positive integers less than m which are relatively prime with m. 

The lower bound in (|4.1|) shows that if n is even, then an n x n circulant sign matrix is much 
more likely to be singular than a sign matrix with independent entries. The first upper bound in 
(|4.2|) implies that the singularity probability is rather lower if n is odd; it is known that 

(4.3) d(n) < n c / loslogn , 

for every n, see [21 Theorem 13.12]. The bound is of course smaller if n has few divisors; for 
example if n = p k for an odd prime p then d{n) = k + 1 = log p n + 1. The second upper bound in 
(|4.2p implies that the singularity probability is extremely small if n has no small prime factors. In 
particular, if n > 3 is prime then 

P[C n is singular] = 2~ n+1 

since this is the probability that all the entries of C n are equal. It would be nice to have a more 
complete description of the dependence of the singularity probability on the prime factorization of 
n. 

Proof. Begin by defining the random polynomial 

n-l 

/(t) = ^Ay, 

j=0 

so that the eigenvalues of C n are /(o»„) for k = 0, . . . , n — 1. Observe first that f(l) and /(— 1) are 
identically distributed and 



(4.4) P[/(l) = ]=P[/(-l) = 0] 



2" n ( n " 2 ) if n is even, 



if n is odd. 

which implies the lower bound in (|4.ip by Stirling's formula. 
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For each k > 0, uj^ is a primitive root of unity of order n/m, where m = gcd(k, n). The minimal 
polynomial of co^ over the rational numbers is thus the cyclotomic polynomial 5> n / m , so /(k>„) = 
if and only if 3>„/ TO is a factor of /. (See e.g. [HI Section V.8] for background on cyclotomic 
polynomials.) Therefore we need only consider the cases when k is a divisor of n, and 

(4.5) P[e n is singular] < P[/(l) = 0] + P = 0] . 

l<fe<n 
fe[n 

A straightforward application of the multidimensional inverse Littlewood-Offord theorem proved 
by Friedland and Sodin [11] yields that 

(4.6) P[/(«£) = 0] <~ 

for some c > when 7^ 0, n/2; the calculations involved in applying the result of [TT] are similar 
to those in the proof of Theorem [TJ The upper bound of (|4.ip and the first upper bound of (14. 2p 
follow by combining (|4.5|) . (14. 4p . and (j4.6H . To bound the number of terms in (|4.5I) in the even 
case it is enough to use the trivial estimate d(n) < 2y/n (divisors of n occur in pairs k,n/k with 
k < \/n) rather than the much more delicate result (|4.3p . 

For n > 3 and d > 2 let -P n)m be the set of polynomials g with rational coefficients of degree at 
most n — 1 such that $ m is a factor of g. With the substitution m = n/k, (|4.5[) becomes 

P[C„ is singular] < P[/ € P n , OT ]. 

Kd<n 
d\n 

Since $ m has degree (p(m), P n ^ m is an (n — y?(m))-dimensional subspace of the (rational) vector 
space of polynomials of degree at most n — 1. Therefore some set of n — <p(m) coefficients of a 
polynomial g suffice to determine whether g £ P n ,mi and so 

P[/€P n , m ] <2~^ m \ 

which proves the second upper bound in (|4.2p . □ 
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